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We develop a Hamiltonian theory for a time dispersive and dissipative (TDD) inhomoge- 
neous medium, as described by a linear response equation respecting causality and power 
dissipation. The canonical Hamiltonian constructed here exactly reproduces the original 
dissipative evolution after integrating out auxiliary fields. In particular, for a dielectric 
medium we obtain a simple formula for the Hamiltonian and closed form expressions for 
the energy density and energy flux involving the auxiliary fields. The developed approach 
also allows to treat a long standing problem of scattering from a lossy non-spherical 
obstacle and, more generally, wave propagation in TDD media. 

1. Introduction 

There is an intrinsic ambiguity in identifying the field energy densities for radiation 
in a time dispersive and dissipative (TDD) medium as described by the linear response 
theory, e.g., in a dielectric medium described by the classical linear Maxwell equations with 
complex valued frequency dependent electric permittivity e(u) and magnetic permeability 
fi(uj). Consequently, there are problems with the interpretation of the energy balance 
equation |T21 §77], |TJ §1.5a], §6.8], [T7j . There were a number of efforts [T7j, [T3] . 
[2H] to construct a consistent macroscopic theory of dielectric media that accounts for 
dispersion and dissipation, based on more fundamental microscopic theories. At first 
sight, it seems that the introduction of a realistic material medium in an explicit form 
similar to ^7j, [T^l, [20] is the only way to model a TDD medium. In fact, that is not so 
and in this paper we describe a consistent macroscopic approach within the linear response 
theory. Full proofs of the statements outlined here will appear in a forthcoming paper [3] . 

A linear response TDD medium is an essentially open dissipative system, which in prin- 
ciple can be obtained by i.) eliminating some degrees of freedom from a more involved 
microscopic theory and ii.) making the approximation of linear response. Stopping short 
of introducing a microscopic theory we ask, is there a conservative extended system which 
exactly reproduces the given linear TDD system after reduction? In |2j we showed that 
indeed such an extension is i.) possible and ii.) essentially uniquely determined, under 
general conditions of causality, power dissipation, and minimality of the extension. Here 
we go further and construct a canonical Hamiltonian for such a conservative extension 
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based only on the given TDD equations — without assumption on the underlying mi- 
crostructure. In particular, we construct such a Hamiltonian for a dielectric medium 
as defined by complex e (u) and [a(uj). The construction given here is not restricted to 
a dielectric medium, however, but holds for TDD systems with a certain mathematical 
structure — eqs. (jHl El HHJ) below — including, in particular, elastic and acoustic media, 
and it can be extended to space dispersive dissipative systems. A somewhat related con- 
struction of the evolution equations for linear absorptive dielectrics was given by Tip [27J . 
The range of validity of the proposed theory is the same as for the linear response, though 
nonlinear generalizations are clearly possible. 

Other important benefits of the approach developed here are: i.) The constructed 
Hamiltonian is an integral of a local energy density, which in the absence of TDD terms 
reduces to the local field energy. This permits us to derive an expression for the energy 
transport for TDD media, ii.) The present formulation allows to treat, in particular, a 
long standing problem of scattering from a lossy nonspherical scatter - analyzed by other 
methods with limited success |18J - by applying the well developed scattering theory 
or conservative systems, see |22|26j and references therein. These applications will be 
discussed in detail in forthcoming work jS], jl]- 



2. Construction of the Hamiltonian 



We consider a system described by two canonical vector coordinates p,q G H, with H a 
real Hilbert space. In the absence of TDD terms, the evolution is supposed to be induced 
by a Hamiltonian A (p, q) of the form 

A (p, q)= l - (K p p , K p p) + X - (K q q , K q q) , (1) 

with closed linear operators K p , K q from H into auxiliary spaces H p , H q respectively. To 
manifest the conservation of energy, it is convenient to consider the evolution of 

f p := K p p G H p , / q := K q q G H q , (2) 

in place of p, q. In the absence of dissipation, these quantities evolve according to 

f P \_f -K\ ff p 



d t |^ P J = I K ^ Q j J , zero dissipation, (3) 
with K := K P K^ a closed linear map from H q to H p . Note that 

^) = ^(ll/pH 2 + ll/ q || 2 ) (4) 

is conserved due to the antisymmetry of the generator in @. 

The electromagnetic field in a non- dispersive inhomogeneous medium may be described 
in this framework, with p = (47r) _1 A (magnetic potential), q = D (electric displacement), 
f p = (2 v /7r)" 1 H (magnetic field), and / q = (2 v /7r)~ 1 E (electric field). Identifying (0) with 
the material relations, we determine the action of the operators K w : 

K p ^(f) = (2 v ^F) / ,- 1 (r) ■ jv x A (f ) J , Kq D{r) = ^-\r) • D(f), (5) 
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where fj,, e are the static permeability and dielectric tensors, assumed real and symmetric. 
We take (4-7r) -1 A and D in the space H = H CUT \ of divergence free vector fields — which 
amounts to a choice of gauge and an assumption of no free charges. To complete the 
picture we define H p , H q to be weighted L 2 spaces with scalar products 

<H,H) Hp = J d 3 fH(r)-/i(r)-H(f), (E , E)^ = J d 3 f E(f) ■ e(f) • E(f). (6) 

As a result 

^ = (2^)Vx, irt = J_p curl) Jf = /i -i(r).Vx, = e-\f)-Vx, (7) 

Zy/TT 

with P cur i the orthogonal projection of (L 2 ) 3 onto H CUT \. 

An alternative formulation of the general system $IJ) is suggested by the example of the 
EM field, namely to consider (j2J) as generalized material relations together with evolution 
equations 

p\ ( ° - K i\ (fi 



In turn, this suggests a natural modification incorporating dispersion and dissipation by 
replacing (J2J) with: 

/■oo 

/ w (i) + / drx w (r)/ w (t-r) = K w w(t), for w = p,q, w = p, q. (9) 
Jo 

The TDD character of @ comes from the operator valued generalized susceptibilities x w , 
w = p, q, the integrals of which explicitly satisfy the causality condition: values of K w w (t) 
depend only on f w (t') for times t' < t. 

Our main result is the following: Assume the susceptibilities x w satisfy the following 
power dissipation condition (PDC), 

Im{Cx w (C)} = ^{CXw(C)-CXw(C) t }>0, w = p,q, for all C = w + i V , V > 0, (10) 

for all £ = u) + ir], rj > 0, where x w is the Fourier- Laplace transform of x w : 
1 f°° 

X W (C) = ^=/ dt^x w (t), w = p,q. (11) 



Then it is possible to construct a Hamiltonian extension to (0p ; which reduces to (QJ) in 
the limit of zero susceptibility, such that the subsystem p,q evolves according to © GJ). 

The extended Hamiltonian A (P, Q) is a function of extended momentum P and co- 
ordinate Q variables, each taking values in a Hilbert space TC D H, and has the same 
structure as ((TJ), i.e. 

A (P, Q) = l - {K V P , /C P P) + l - (/C q Q , /C q Q) , (12) 
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with /C p , /C q closed operators from Ji to H p D H p , Tt q D i? q , which extend K p and -K" q 
respectively (see f)48|l below). 

Before presenting the general construction, let us illustrate it with the example of a 
linear TDD dielectric medium, described by the macroscopic Maxwell equations without 
external charges and currents 

dtD = V x H, SfB = —V x E, V ■ B = 0, V ■ D = 0, (13) 
in units with c,e ,fi = 1. Here 

D = E + 4vrP, B = H + 47rM, (14) 



with the polarization P and magnetization M given by linear response 

poo poo 

P(r,t)= drx E (r,T)E(r,t-r), M(f,t) = dr x H (r, r) H {r, t - r) . (15) 



The electric and magnetic susceptibilities should satisfy the PDC (jlUj) for each r, and for 
simplicity we take them to be real valued scalars. (The frequency domain susceptibilities 
may nonetheless be complex.) 
Motivated by |2| and the Lamb model (see Fig. ^ below), we introduce canonical vari- 
ables 



P = ((4^-^(0, 9 E (r, s), y> H (r, s)) , Q = (D(f), <p E (f, s), H (r, s)) 



(16) 



with A, D G H cm \ and auxiliary vector fields ip F , 6 F , F = E, H, which are functions of f 
and an auxiliary coordinate — oo < s < oo. For these variables we define a Hamiltonian 



A(P,Q) = T(P)+U(Q) 



with, 



i r _ A(r) r°° 

T(P) = 2j d¥ 2 v^V x -A-2 - J ds<; H (f,s)<p n (r,s) 

"OO 



+ 



\ J d 3 rj™ds [\0 E (r, S )\ 2 + \d s cp H (r,s)\ 2 ] 



(17) 



;is) 



U{Q) = -J d 3 f D(f) - y ds <j E (r, s)cp E (f, s) 

d S [|^ E (f, S )| 2 + |0 H (r>)| 2 ] 



i /dV / x 



(19) 



where q F , F = E, H, are scalar functions to be specified below. 
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The resulting Hamilton equations of motion for the extended Maxwell system are: 

dt^- = ~ [D(r, t) - 2^ fas, Ve) s (r, t)] , (20) 
d t<PH ( r > s ) = -On (r, s) , 

d t 6 E (r, s ) = J- ?E (r, s) [D(r, t) - 2^ ( ?B , Ve ) s (r, t)] + «9 s Ve (r, s) , 

Zy/TT 

d t D(r, t) = V x [V x A(r, t) - 2^ (? H , <Ph> s (r, *)] , (21) 
d t H (r, s) = ^ H (r, s) [V x A(r, t) - 2^ (? H , ¥>h>* <)] - ^Vh (r, s) • 

^t^E ( r > s ) = 0e (r, s) , 
where 

(•,•>. = / •• da. (22) 
We make the identifications: 

B(r,t) = VxA(r,t) (23) 
H(r,t) = B(r,*)-2V^H,V> H ) S M) = B (r, t) — AttM. (r, t) , (24) 
E(r,t) = D(r,t)-2v^(^ E ,^ E ) s (r,t) = D (r, £) - 4ttP (r, t) , 

with 

<? H , ¥>h> s (r, t) = 2^ M (r, *) , ( fE , c^ E ) s (r, t) = 2^ P (r, t) , (25) 
resulting in the following equivalent system of extended Maxwell equations 
d t H (r, t) = — V x E (r, t) + 2^ ( ?H , 6> H ) S (r, t) , (26) 
^t^H ( r , s ) = (r, s) , d t E (r, s) = ^ E (r, s) E (r, t) + <9 S V E (r, s) , 

2V7T 

<9 4 E (r, t) = V x H (r, *) - 2^ ( ?E , E ) s (r, t) , (27) 

<9 t H (r, s) = -77^ (r, s) H (r, t) - d 2 s <p u (r, s) , d t ip E (r, s) = E (r, s) . 
2V 71 " 

Combining the first order Hamilton equations of motion for ip F and 6 E into a single 
second order equation for cp F , F = E, H, we obtain a driven wave equation 

{#-#}M^M) = 2^?F(* , ,a)F(r' ) f) > F = E,H. (28) 

Assuming y> F to be at rest (<p F = d t <p F = 0) in the distant past, the solution to (J28|) is 
given by 

poo rs+r 

<p F {r,s,t) = ^j= I dr I da^ F (f,a)F(f,t-r), F = E, H, (29) 



o 
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implying with (|T5|) and the following expression for the susceptibilities 



^ pco rs+t 

X F (r,t) = — ds d<7<r F (f»<r F (f», F = E, H. (30) 

O 71 " J -00 Js-t 

The key fact is that, due to the power dissipation condition (|10|). it is possible to invert 
(f3~0|) and write <j F as a function of x F . An explicit solution is 

? F (f, *) = -==/ dcue^V^ Im XF (r> + i0), F = E,H. (31) 

V27T J-oo 

Note that <r F is real and invariant under sh>-s. 

The above discussion implies the following result on the extended system: Let the 
Hamiltonian \1T[ [75|. 6e groen with q-p defined by for Xf which obey Then 
for any solution to the Hamilton equations of motion with if F ,8p — > as t — > —00, the 
variables D(r,i) and B(r, t) = V x A(r, t) evolve according to the macroscopic Maxwell 
equations 173103)) . 

Based on the constructed TDD Hamiltonian ()17jl . we obtain an expression for the 
energy density of the EM field and the medium 

S (r,t) = i- {|E| 2 + |H| 2 } (f,t) + \ {\\d.<pX + ||0h||' + II^eII' + ||*e|£} (r,t) , (32) 
where 

a 1 '-' 2 ds. (33) 
This results in the conservation law 

d t £ + V • S = 0, (34) 
with the familiar Poynting vector for the energy flux 

S(f,t) = i-E(f,t)xH(f,t). (35) 

These identities follow from (jl7H19)) and the general theory of Hamiltonian fields 

When the interaction ^ is set to zero, the EM and auxiliary fields decouple and (J32|) 
reduces to 

S (f, t) = Sem (r, t) + S s (f, t) (36) 
with the energy density of the EM field 

£em (r, t) = |D| 2 (f, t) + J- I V x A| 2 (f, t) , (37) 
and the energy density of the auxiliary fields 

£s (r,t) = \ {H^hII 2 + ||0 H |L + II^eIIJ + II^eII:} (r,t) . (38) 
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Subtracting (|55|) from (|3~2"j) gives the energy shift due to the interaction of the EM field 
and the matter 

68 (f,t) = -L {|E(r,t)| 2 - |D (r,t)| 2 + |H (r,t)| 2 - |B (r, t)| 2 } . (39) 

671 

In general, it is not possible to give an expression for the energy density £$ (r , t) of the 
medium in terms of the instantaneous EM fields E(r, t) and H(r,t). However, using (J291 
l3*U|) and the equations of motion we have calculated that 

dt \ {ll^pll 2 + ||0 F || 2 } (r,t) = [d t (^cp F )(f,t)] • F(f,t) (40) 

= [d t P(f,t)}-E(r,t), F = E, (41) 

with a similar expression for F = H. 

The result (14 lj) is the usual expression for the rate of change the density of EM energy 
stored in a dielectric. For a wave packet E(f, t) = Re {e~ laj °*Eo(r, t)} with Eo(r, t) a 
slowly varying function of t, we have 

d t P(f,t) « v^Re |-L;oXE(^^o)e- i " of Eo(F,t) + ^"^'^^ e-^^Eo^, t)| , (42) 
and thus 

[d t P(f,*)]-E(r,t) 

« ^Re |-icoX E (r^o)|Eo(f,t)| 2 + d ^^° } [frEo(r,t)] ■ E*(r,o| 

+ terms with a factor of e~ itJ °* or e^ '. (43) 

If we consider the time averaged power density, averaged over a time scale longer than 
1/coo, the terms with oscillatory factors are very small and we have (with denoting time 
averaging) : 



[d t I>(r,t)]-E(r,t) 
fn { 1 



u; Rex E ( r ? ^oJ 



d^ 



^|E (f,t)| 2 + ^oImx E (r,^o)|Eo(f,t)| 2 ^ , (44) 



where we have assumed for simplicity that the slowly varying function E (r*, t) is real. 
The first term on the r.h.s. of (pij) is a total derivative, which when integrated gives the 
Brillouin formula for the time averaged energy density in a material medium (see [T2"l 
Section 80]): 
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|E (f,t)| 2 , (45) 



which neglects losses entirely, and hence is useful only if Imx E (r, uq) is zero or so small as 
to be irrelevant. The second term of (j44j) is strictly positive and gives the contribution 



2 



cj Imx E (r> ) / |E (f,r)| 2 dr ( 1.6) 
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to the time averaged material energy density, which is non-decreasing in time and incor- 
porates losses in an approximate way. 

For a general system of the form (|5| I9| ITU|). a TDD Hamiltonian can be constructed in 
the same way. We define canonical variables 

P = (p, 9 q , <p v ), Q= (q, w 6 P ) , (47) 

with 9 w (s), <p w (s) functions of an auxiliary coordinate — oo < s < oo, taking values in the 
Hilbert spaces H w , w = p, q. If Xw{t) = Xw(*) f° r an ^ > 0, the Hamiltonian is of the 
form (|12p. with 

/C P P = \K p p- J dsg p (s)(p p {s), 6 V d a (p p J , (48) 
where 

e„H = (27r)- 1 / 4 v ^m^>TT0), w = p, q, (49) 

with yfi the operator square root. 

In particular, we can handle it this way: i.) Non-isotropic media, provided the tensors 
X w are real symmetric. (Gyrotropy could in principle be handled by a more involved 
construction with terms mixing momenta and coordinates.) ii.) Space dispersion, in 
which case terms depending on V(£> and V0 appear in the Hamiltonian. Details of the 
abstract construction and further examples will be given in forthcoming work [Hj. 

3. Discussion and comparison with prior work 

The need for a Hamiltonian description of a dissipative system has long been known, 
having been emphasized by Morse and Feshbach [T^ Ch. 3.2.] forty years ago. They 
constructed, for a damped oscillator, an artificial Hamiltonian based on a "mirror-image" 
trick, incorporating a second oscillator with negative friction. The resulting Hamiltonian 
is quite un-physical: it is unbounded from below and under time reversal the oscillator 
is transformed into its "mirror- image." The artificial nature of this construction was 
described in j!9| Ch. 3.2.]: "By this arbitrary trick we are able to handle dissipative 
systems as though they were conservative. This is not very satisfactory if an alternate 
method of solution is known..." 

The Hamiltonian we construct for TDD media can be viewed as a quite general "satis- 
factory solution" to the problem posed in |THJ Ch. 3.2.] since we do not introduce negative 
friction and, in particular, we do not make use of "mirror-images." Instead we couple a 
given TDD system to an effective model for the normal modes of the underlying medium. 
For the combined system we give a non-negative Hamiltonian with a transparent inter- 
pretation as the system energy. As regards the underlying microscopic theory, this is an 
effective Hamiltonian for those modes well approximated by linear response. 

The evolution equations of the proposed theory come from a Hamiltonian and are thus 
time reversible. Nonetheless, an irreversible motion of the TDD system stems from the 
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Figure 1. The Lamb model, introduced in ^0] to describe radiation damping, is a point 
mass attached to an infinite elastic string and a Hook's law spring. The point mass evolves 
SIS db classical linearly damped oscillator. 



infinite heat capacity of the auxiliary system. This is demonstrated in its simplest form 
by the damped harmonic oscillator 

md 2 t q(t) + jd t q{t) + kq(t) = 0, (50) 
which results from the TDD system 



9 t [ p \ = [ ; q ) , r pw 1 uo f / + r \ = (51) 



k\ (IS) +7/ °°/p(*-T)dr\ f^P(t)\ 



qj \f P J ' V /q(*) J \Vkq(t)J 

with m, fc,7 > 0. The construction presented above reproduces in this simple case a 
model due to Lamb in 1900 10 - see Fig. ^ — in which the energy of an oscillator 
escapes to infinity along an attached flexible string. The theory proposed here illustrates 
that, from the standpoint of thermodynamics, dissipation in classical linear response is 
an idealization which assumes infinite heat capacity of (hidden) degrees of freedom. 

In general, the auxiliary system, described by the fields {y2 p , (p q , 9 p , 9 q }, is governed by 
a Hamiltonian B.^((p, 9) of the following simple and universal form 

1 r°° r i 

H hb (^) = -[ \\e(s)\\l v9 ^ + \\d.<p{8)\\ 2 Hit9Ha ds. (52) 



The Hamiltonian Hhb(y?, #) in ()52j) is a canonical heat bath (justifying the index hb) as 
described in [HJ Section 2], Section 2]. 

The physical concept of an ideal or canonical heat bath originates in thermodynam- 
ics. Statistical mechanical models at the mathematical level of rigor were introduced, 
motivated, and described rather recently (to our best knowledge), see [3 Section 1], 
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Section 2] , [221 Section 2] and references therein. The consensus among the references, on 
the basis of statistical mechanics, is that the generator of motion for a canonical heat bath 
must be ix a self-adjoint operator with absolutely continuous spectrum with no gaps, i.e. 
the spectrum must be the entire real line K, and the spectrum must be of a uniform 
multiplicity. These requirements lead to a system unitarily equivalent to the universal 
form Hamiltonian ([52)1 . [HI Section 2]. 

General statistical mechanics considerations indicate that for a system to behave ac- 
cording the thermodynamics it must be properly coupled to the heat bath. In particular, 
the coupling should involve all modes of the heat bath, |7J Section 1], Section 2], [213 
Section 2]. In our Hamiltonian setting the coupling is essentially (K w w,T w ip w ), which is 
the dipole approximation, [23 Section 1,2], and the condition of coupling to all modes is 
the following constraint on the dissipation: 

[J Ran Imx w O) = H w , w = p, q. (53) 




There is some relation in spirit between our theory and a recently proposed hydro- 
dynamic theory (HT) ^H[> 111 Both theories are self contained, macroscopic, and 
make no assumption on the underlying microstructure. However the theory proposed 
here, unlike the HT, makes no use of parameters other than the susceptibilities of linear 
response theory. Furthermore, the present theory is truly conservative, with dissipative 
effects modeled by effectively irreversible energy transport to an auxiliary system, which 
may be conceived of as constructed from flexible strings. In contrast, the HT makes use 
of explicitly dissipative, nonconservative equations similar to those of Navier-Stokes. 

A deeper relation can be found between the approach described here and the well known 
dilation theory which, in certain cases, provides a treatment for dissipation and resonance 
phenomena. We first give a brief account of the dilation theory, based on B. Pavlov's 
extensive review [2HJ as well as his more recent work [21]. For more detailed exposition 
on the subject we refer the reader to [22], [21] and references therein. 

The dilation theory was the first rather general approach to the construction of a 
spectral theory for dissipative operators. It is based on an abstract version of the Lax- 
Phillips scattering theory fj] which assumes that there are i.) a dynamical unitary 
evolution group U t = e int in a Hilbert space H where Q is a self-adjoint operator in H; ii.) 
an "incoming" subspace D_ C H invariant with respect to the semi-group U t , t < 0, and 
an "outgoing" subspace D + C H invariant with respect to the semi-group U t , t > 0. The 
invariant subspaces (also called scattering channels) D± are assumed to be orthogonal. 
Then one introduces the "observation" subspace K = HQ(D_ © D + ) which is coinvariant 
in the sense that the restriction of Ut, t > 0, to K is a well defined semigroup on its own, 
namely for t > 

Z t = Pi<U t \ K = e iBt , where P K is the orthogonal projection on K. (54) 

In many interesting cases the generator B of the semigroup Z t is dissipative, i.e. Im B > 
or, even, ImB > 0. So the relation (|54j) provides an interesting scenario within the Lax- 
Phillips scattering theory for the rise of a dissipative operator. The dilation theory yields 
the spectral theory through the construction of generalized eigenmodes (the scattering 
theory) provided, of course, the conditions discussed above are satisfied. 
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Looking at the dilation theory from the point of view of open (dispersive and dissipative) 
systems, one can ask if the theory allows to find the unitary group U t = e lUt or, equiv- 
alently, the self-adjoint operator Q being given the dissipative operator Bl The answer 
is positive for rather large class of dissipative operators B. For example, if B = Q + ia, 
where f2 is self-adjoint and a > is bounded, a unique minimal dilation and its eigen- 
modes can effectively constructed, [2U Theorem 3]. 

Hence when the dilation theory applies it provides a solid foundation for spectral studies. 
Unfortunately, the dilation theory does not apply to many important physical problems 
simply because its initial assumptions on the nature of the dissipation are too restrictive. 
For systems described by evolution equations ©-© the dissipation always comes with 
the dispersion, and the dilation theory does not apply. Indeed, the most general form for 
a linear causal time-homogeneous open system, as analyzed in 0, is 



where H Q is a Hilbert space, m > and Q are self-adjoint operators in H , f (t) is an 
external force and a (t) is the friction function, subject to a power dissipation condition 



Only in the very special case of \55)) when the friction is instantaneous (Markovian) , i.e. 
a(t) = ao5(t), can one use the dilation theory as in \2J\ Theorem 3]. For many well 
studied dielectric media, such as Lorentz or Debye, not to mention media with generic 
frequency dependent electric susceptibilities, the relevant friction functions are not in- 
stantaneous. For such systems one must use a more general approach, such as developed 
in [2] and extended here and in [3]. 

It is interesting to point out, however, that coupling to a canonical heat bath (jo^j) may 
be interpreted, as in the Lamb model, as "attaching an elastic string" at any point of loss. 
The attached strings are analogous to the scattering channels of the dilation theory, with 
the only difference being that our "strings" are coupled in more general ways than in the 
dilation theory. Thus, one can view our TDD Hamiltonian as a natural generalization of 
the constructions of the dilation theory. 
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